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ABSTRACT 

We report the first detection of flux variability in the most luminous X-ray source in 
the southern ring of the Cartwheel galaxy. XMM-Newton data show that the lumi- 
nosity has varied over a timescale of six months from £o.5-iokcV ~ 1-3 x 10 41 ergs -1 , 
consistent with the previous Chandra observation, to L u .5-iokcV ~ 6.4 x 10 40 ergs -1 . 
This fact provides the first evidence that the source is compact in nature and is not 
a collection of individual fainter sources, such as supernova remnants. The source has 
been repeatedly observed at the very high luminosity level of Lo.5-iokcV ~ 1-3 x 10 41 
ergs -1 for a period of at least 4 years before dimming at the current level. It rep- 
resents then the first example of an accreting object revealed in a long lived state of 
extremely high luminosity. 
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1 INTRODUCTION 

Very luminous off-nuclear X-ray sources were discovered 
, in nearby galaxies with the Einstein satellite (Fabbiano et 
■ al. 1989). They were named Ultra-Luminous X-ray sources 
(ULXs) based on their X-ray luminosities, much higher than 
the Eddington limit for a solar mass black hole (Lx ~ 
1.4 x 10 38 ergs -1 ). These luminosities may reflect beamed 
emission from an accreting stellar mass compact object, or 
super- Eddington emission, or isotropic accretion onto an in- 
termediate mass black hole. The brightest objects, those 
with Lx ^ 10 41 ergs -1 , sometimes termed Hyperluminous 
X-ray sources (HLXs; see Matsumoto et al. 2001; Kaaret et 
al. 2001), are even more intriguing, since their luminosities 
are closer to that of active galactic nuclei, requiring a bigger 
engine, stronger beaming or even extreme super-Eddington 
regimes. The issue of the physical interpretation of ULXs 
and HLXs is still quite open. 

An extraordinary example of HLX is the source N.10 
detected in the narrow, gas-rich star-forming ring of the 
Cartwheel galaxy with isotropic luminosity of Lo.5-iokcV ~ 
1.3 x 10 41 ergs -1 (Wolter et al. 1999; Wolter & Trinchieri 
2004 - hereafter WT04; Gao et al. 2003). This is the brightest 
of a number of individual sources that also appear to reside 
in the ring, all classified as ULXs, based on their isotropic lu- 
minosities in excess of Lo.s-iokov = 3xl0 39 ergs -1 (WT04). 
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The spatial association of the N.10 HLX and the ring- 
ULXs with HII complexes and young star-forming clus- 
ters suggests a physical link, thus providing an invaluable 
in-depth probe of the young stellar population currently 
present in the Cartwheel outer ring. The physical nature 
of the ULXs in the Cartwheel is not clear yet: they could be 
genuine single sources, hence accreting compact objects, or 
unresolved collections of supernova remnants. A way to dis- 
entangle this puzzle is to look at their variability. Variability 
in ULXs, covering timescales of months to a few years, has 
been often encountered (e.g., in the Antennae; Fabbiano et 
al. 2003) and was taken as evidence that they represent an 
apparently bright state of accreting high mass X-ray bina- 
ries (HMXBs) caught in some peculiar evolutionary stage 
(King 2002). Young, short lived, HMXBs hosting a neutron 
star or a stellar-mass black hole are common in star-burst 
galaxies and their number is likely to be linked with the 
star-formation rate of their hosting galaxy (Fabbiano 2005). 

It has been proposed that HMXBs share a universal X- 
ray luminosity function that extends up to luminosities of 
10 40 ergs -1 , characteristic of ULXs (Grimm, Gilfanov, & 
Sunyaev 2003). WT04 have shown that the X-ray sources 
seen in the Cartwheel follow the same luminosity function 
and derive a star formation rate of 20 M©yr -1 in agreement 
with recent radio and Far Infrared estimates (Mayya et al. 
2005). The only outlier appears to be source N.10, brighter 
than expected at the bright end of the X-ray Luminosity 
Function (XLF) by a factor of 3. WT04 noticed that a higher 
cutoff would account for this excess in the Cartwheel XLF, 
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which would otherwise highlight a different physical origin 
for this exceptional HLX. 

At present, the only confirmed HLX reported in the lit- 
erature is the source X-l in the starburst galaxy M82, with 
isotropic luminosity of a few x 10 40 erg s -1 (Ptak & Griffiths 
1999: ASCA baseline flux; Strohmayer & Mushotzky 2003; 
Dewangan et al. 2006: XMM-Newton observations). X-l in 
M82 showed a burst, up to 10 41 ergs -1 , lasting about a 
month during ASCA observations (Ptak & Griffiths 1999) 
and again with Chandra (Matsumoto et al. 2001; Kaaret 
et al. 2001). Spectral and temporal variability points to a 
2OOM0 intermediate mass black hole as the accreting ob- 
ject (Dewangan et al. 2006). If confirmed by additional ob- 
servations, X-l is the prototype of a "new" class of black 
holes, bridging the gap between those of stellar origin and 
the very massive black holes hosted in galactic nuclei (Miller 
& Colbert 2004). Other candidate HLXs are those in NGC 
7714 (L x = 6. X 10 40 erg s -1 in the high state; Soria & 
Motch 2004; Smith, Struck & Nowak, 2005), in NGC2276 
{Lx = 1.3 x 10 41 erg s -1 in the 0.5-10 keV band; Davis & 
Mushotzky 2004) and in MCG-03-34-63 (ULX 1 with peak 
luminosity of Lx = 1.3 X 10 41 erg s -1 in the 0.5-7 keV band; 
Miniutti et al. 2006). The majority of these objects have 
been observed with XMM-Newton only, and with limited 
time coverage; therefore neither a secure physical explana- 
tion nor a clear association with the host galaxy can be put 
forward. 

The Cartwheel source N. 10 is possibly the brightest 
HLX known but the lack of detailed information on its vari- 
ability has prevented the exclusion of an extended nature. 
The new data that we present here have now confirmed its 
compact nature. 

The X-ray image of the Cartwheel shows, besides a 
number of point sources, the presence of hot (fe_eT=0.2 
keV) gas in the ring with a luminosity of Lg a ^_ 2kcV — 
3 x 10 40 ergs -1 (WT04) and suggests that more gas might 
be present in the near environment of the galaxy. To better 
study all of these components we have obtained new XMM- 
Newton observations. Here we concentrate on the proper- 
ties of the HLX, leaving a more comprehensive presentation 
of the galaxy as a whole to a forthcoming paper. For the 
Cartwheel we use a distance D=122 Mpc (calculated from 
Amram et al. 1998, with Ho = 75 km s -1 Mpc -1 ) corre- 
sponding to a linear scale of 0.592 kpc/arcsec. 



2 XMM DATA 

The Cartwheel was observed with XMM-Newton on the 14- 
15th of December, 2004 (36 ksec; dataset [101]) and on the 
21st-22nd May, 2005 (60 ksec; dataset [201]), with the pn 
and MOS instruments operating in full-frame mode with the 
thin filter applied. Standard rejection criteria were applied 
to eliminate high background periods that are fortunately 
relatively short. The net exposure times, after data clean- 
ing, reduce to 29(24) and 50(42) ks for MOS(pn), a ~ 20% 
reduction in time. 

The two datasets were analyzed independently. We have 
used the standard XMM Science Analysis System tasks to 
prepare the data and produce images and spectra. Our pri- 
mary interest, source N.10 in the Chandra list (WT04), falls 
over a bad column in MOS2 in the first observation, and 



in MOS1 in the second one. We will therefore consider only 
MOS1 and pn for the [101] observation and MOS2 and pn 
for the [201] observation. 

We have checked consistency between the two obser- 
vations by using different regions within and outside the 
Cartwheel. The aspect solutions of the two dataset are very 
similar, so that source positions do not differ between [101] 
and [201]. Bright field sources are expected to be mostly 
AGN, therefore flux variations of the order of a factor of 2 
are common, so no single source can be used as a reference 
point. However, we have confirmed that overall the count 
rates are constant for a large number of sources between the 
two observations. 

2.1 Image 

We compare the pn images from the [101] and [201] observa- 
tions, smoothed with an adaptive Gaussian kernel (package 
csmooth from Ciao 3.3). In Fig.^we plot the two images side 
by side. We overplot for reference the positions of Chandra 
sources, that are not at the peak of the XMM-Newton posi- 
tions. In fact, no correction for the different aspect solutions 
of the two satellites has been attempted, however distances 
between Chandra and XMM-Newton peaks are below as- 
pect uncertainties . Although the two images should not be 
used for a quantitative comparison since the two observa- 
tions have different lengths and should be properly normal- 
ized and corrected for possible background differences, the 
graphical comparison shows variability in different areas. In 
particular, it is evident that the HLX is no longer the bright- 
est source in the second observation. From the comparison 
between the count rates in the neighbouring sources, we have 
determined that the source next to it, which corresponds to 
Chandra sources N.13 & N.14, has not varied between the 
two observations (the count rate in the second observation is 
at most 15% higher than in the first one). If we assume that 
this is constant, then the source to the NW, corresponding 
to N.16 & N.17, is also constant (the same 15% increase in 
the count rate), but the SE source (N.7 & N.9) has faded to 
about half its strength in the second observation. Note that 
the darker colours in the image simply reflect the higher 
statistics available in the second exposure, which is about 
twice in length. 

2.2 Spectrum 

We extract spectra from a region of radius 10" cen- 
tred about the peak of N.10 in [101] and background from 
a nearby circular region devoid of sources. Appropriate re- 
sponse matrices for spectral analysis where generated using 
the SAS tasks arfgen and rmfgen. To improve on the statis- 
tics, we have binned the data so that each bin has a sig- 
nificance of at least 2a. We report in Tabled the total net 
counts and exposure times in seconds, after cleaning, for this 
extraction region. 

The extraction radii for XMM-Newton are smaller than 
customary, due to the presence of many surrounding sources 
in the Cartwheel ring. Nonetheless they include a portion of 
the ring, so we expect a fraction of the diffuse underlying 
gas component to contaminate the HLX spectrum. 

We first fit the [101] spectrum with a simple model, 
i.e. an absorbed power law. This results in an acceptable 
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Figure 1. A smoothed grey scale image of pn data for the southwestern ring from the first observation [101], left, and the second [201]. 
right. The positions of the brightest Chandra sources are indicated. No correction for different aspect solutions has been attempted. 
Distances between Chandra and XMM-Newton peaks are below aspect uncertainties. In the [101] data, N.10 is clearly the brightest 
source; it is unresolved by XMM-Newton. In the [201] observation, which is longer, more details and more sources appear. The detection 
of a faint source at the position of N. 10 is difficult due to the close vicinity of the brighter sources N. 13&14. Notice also the relative 
intensity of sources N.7&N.9 and sources N.16&N.17 which is varied between the two observations. Sources N.7&N.9, N.13&N.14, and 
N.16&N.17 cannot be separated by XMM-Newton. 



Table 1. Log of XMM-Newton Observations 



Name 


Date 


Instrument 


Net Counts 11 


Exp. Time 








(0.5-10 keV) 


sec 


N.10. 


12/14/2004 


MOS1 


80.8±9.5 


29,583 




[101] 


pn 


202.1±14.6 


24,418 


N.10. 


05/21/2005 


MOS2 


87.2±9.9 


49,669 




[201] 


pn 


237.9±15.9 


42,277 



in a 10" radius, centred on RA(2000) = 00 h 37 m 39*38 and Dec(2000) 



-33° 43' 23'.'08, see text. 



fit (x 2 = 15.7 for 23 d.o.f) with a slope r = 1.75 ± 0.25, 
low energy absorption due to an intervening column with 



N H = 1.4 x 10 21 cm" 2 and flux f . 5 



io = 6.9 x 10" 14 erg 
. This is consistent with the Chandra results for 



cm s 

source N.10 suggesting that the HLX is the dominant source 
of emission. However, due to the larger area considered and 
larger PSF, we expect some contribution from the under- 
lying ring emission. As will become clearer later, this as- 
sumption is also relevant for the comparison between the 
two XMM-Newton observations. We have therefore consid- 
ered a spectrum that includes all three components derived 
in the Chandra data: 1) gas, 2) unresolved binary sources, 3) 
HLX (from WT04). 

Given the limited quality of the data, we have fixed the 
parameters of components 1) and 2) to the same values of 
WT04. For the relative normalizations we used the Chandra 
values, properly rescaled to the area covered by the present 
region (1/10 of the emission detected by Chandra in the 



Table 2. Fluxes of different fit components 



Component 


Flux 


Flux 




(0.5-2 keV) 


(2-10 keV) 


HLX [101] 


2.4 X 10 -14 


5.0 x 10 -14 


HLX [201] 


1.2 x 10~ 14 


2.3 x 10~ 14 


gas 


2 x IO" 15 




unresolved binaries 


1 x 10~ 15 


1.3 x 10~ 15 



whole ring). The HLX normalization is left free, while F is 
fixed at the Chandra value. The details of the model are 
shown graphically in the rightmost panel in Fig. 2. 

The HLX is indeed the brightest component in the 10" 
radius region; see Table El for fluxes of the different compo- 
nents. The intrinsic luminosity of the HLX is Lo.5-2kcV — 



4.6 x 10*" ergs -1 and L 2 - 



.7 x 10 4U ergs" 
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channel energy (keV) channel energy (keV) channel energy (kev) 

Figure 2. Left: spectrum of [101], fit as described in the text. Centre: spectrum of [201] with fit from [101]. Right: Unfolded spectrum, 
from pn data only for clarity, showing the three components at the best fit values. The difference between top [101] and bottom [201] is 
only the normalization of the 'HLX' component [and the binning scheme] . 



We then used the same extraction region, and the same 
binning scheme for the [201] dataset. We obtained about the 
same number of net counts in spite of the longer observing 
time. Again, the single power law gives an acceptable fit 
(X 2 = 18.8 for 21 d.o.f.) with parameters consistent to [101], 
and a flux about half. If we apply the same complex spectral 
model of [101], with the same parameters, we find a large 
discrepancy with the data (x 2 > 170 for 21 d.o.f.) as shown 
in Fig. (Centre). In particular we notice that the points 
above 0.8 keV are systematically lower than the model. If 
we make the reasonable assumption that both the diffuse 
gas and the unresolved point source components have not 
varied, and let only the normalization for the HLX compo- 
nent vary, we obtain a good fit for a normalization that is 
about 1/2 that of [101] (x 2 = 21.13 for 23 d.o.f.). Given the 
data quality and the complexity of the spectral model, test- 
ing a spectral variation in this component is unrealistic. We 
therefore cannot comment on spectral variations between a 
high and a low state expected for binary sources. However 
we can reasonably state that the HLX has dimmed by at 
least a factor of two in the 6 months between the two obser- 
vations. The unabsorbed flux of N.10 in this second observa- 
tions at the formal best fit values is reported in Tabled The 
corresponding luminosity is at most Lo.5-2k c v = 2.0 x 10 40 
- 1 and i 2 -iokov = 3.8 x 10 40 



ergs 



2.3 Light curve 



erg s 



The evidence for variability during the XMM-Newton 
observation prompted us to look at the long term behavior 
of the source. We construct the light curve of source N.10 
by using all the available datasets (Fig. [3J • We extract net 
counts from a 10" radius region in the 0.3-2.5 keV band, to 
match the limited energy band of the HRI, and use a nearby 
large circle devoid of sources for background. We convert 
count rates to flux by using the Chandra fit with a power law 
and conversion factor from PIMMS (see WT04). Of course 
cross calibration uncertainties are possible between different 
instruments, and we cannot be sure that the shape of the 
emission was the same at all times, however these uncertain- 
ties are probably small in the band considered. 



HRI 

i i i i | i i i 



CHANDRA 
I I I I | I I I 



XMM 
— I — I — h- 



_ i i i i 



4.9X10- 1 5xl0 4 5.1xl0 4 5.2x10* 5.3x10" 5.4xl0 4 
Time (JD - 2400000) 



Figure 3. Long term light curve in the soft (lower panel) and 
hard (upper panel) energy band, over an interval of about 10 
years. The two XMM-Newton points that define the variation 
are not subject to cross-calibration uncertainties. All fluxes are 
computed with the same spectral shape, i.e. power law with T = 
1.6 and N H = 3.6 X 10 21 cm -2 ; see fit to Chandra data in WT04. 
Right axis reports luminosities computed assuming the Cartwheel 
distance. 



Dates, instruments, count rates and fluxes in both the 
soft and hard band, extrapolated from the same model, are 
reported in Tab. |H] From inspection of the light curve we 
deduce that the HLX was in a brighter state from 2001 to 
2004. The HRI errorbar is consistent at the lowest level with 
a "ring-only" contribution, but an enhancement of the emis- 
sion was visible in the area in the HRI data, to indicate that 
the source was probably also in this bright state since 1994 
(see Wolter et al. 1999). The source became significantly 
fainter in the six months between the 2004 and 2005 XMM- 
Newton observation, assuming a "ULX" status. 
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Table 3. Summary of all observations of N.10 



Date 


JD 


Instrument 


Count rate 


Flux 


Flux 








(0.3-2.5 keV) 


(0.5-2 keV) 


(2-10 keV) 



1994 Dec 9-23 2449696-2449710 ROSAT HRI 3.79 ± 1.08 x 10" 4 2.6 x 10~ 14 5.6 x 10~ 14 

2001 May 26-27 2452056-2452057 Chandra ACIS-S 5.30 ± 0.26 x 10" 3 3.0 x 10~ 14 6.7 x 10~ 14 

2004 Dec 14-15 2453354-2453355 XMM-Newton pn 7.8 ± 0.57 x 10~ 3 3.4 x 10~ 14 7.2 x 10~ 14 

2005 May 21-22 2453512-2453513 XMM-Newton pn 4.2 ± 0.33 X 10~ 3 1.8 X 10~ 14 3.9 X 10~ 14 



3 DISCUSSION 

Spectral properties and variability in some ULXs (e.g. Mak- 
ishima et al. 2000) suggest that we are witnessing accretion 
onto a compact object, in a binary system. However, no uni- 
versal model exists. ULXs might be stellar mass black holes 
with anisotropic X-ray emission due to mechanical beam- 
ing (King et al. 2001), or relativistic beaming from a jet 
(Mirabel & Rodriguez 1999; Kording et al. 2002), or stellar 
mass black holes accreting at super-Eddington rates (Begel- 
man 2002). The most challenging model for the HLXs, given 
their extreme luminosities, is that of a binary system hosting 
a 10 2 ~ 4 M black hole (e.g. Colbert & Mushotzky 1999). In- 
termediate mass black holes may form from the collapse of 
very massive stars born through stellar runaway collisions 
in dense star clusters (Portegies Zwart & McMillan 2002; 
Gurkan et al. 2004). In such a young environment, the in- 
termediate mass black hole may gain a massive donor star 
( > 20 M Q ) through tidal events (Baumgardt et al. 2005) 
or dynamical interactions (Blecha et al. 2006). This system 
will be able to sustain luminosities as high as io 40-41 ergs -1 
as seen in the very bright ULXs (Patruno et al. 2005; Mad- 
husudhan et al. 2006). Spatial association of X-ray sources 
with young star cluster has been searched for a number of 
ULXs (Kaaret et al., 2004) but has not been firmly estab- 
lished yet. The whole ring of the Cartwheel consists of bub- 
bles and condensations (Struck et al. 1996) and the neigh- 
bourhood of N.10 is no exception. Given the distance of the 
galaxy any small misalignment between X-ray and optical 
(HST) positions implies kpc scale distances, so a precise de- 
termination of the optical counterpart is hard, without a 
proper absolute cross-calibration of the two images. In any 
case, the association with an environment of massive and 
young stars is almost certain. 

The XMM-Newton observatory has revealed a factor of 
two dimming in the flux of source N.10 in the Cartwheel, 
over a timescale of ~ 6 months. Although this kind of vari- 
ability has been observed in other ULXs, in the N.10 case it 
provides the first evidence against the hypothesis that the 
source is a chance superposition of fainter sources such as 
young supernova remnants, and suggests that it is a truly 
compact source. The estimated age of the southern ring 
(< 10 Myr), and the lack of radial spread of its sources, 
indicate that source N.10 is closely linked to the active star 
forming episode and its youth suggests that the hypotheti- 
cal intermediate mass black hole present in N.10 had a high 
chance of capturing the massive companion star through a 
tidal event (Baumgardt et al. 2005) . The decline in flux by a 
factor of two brings the source down to the luminosity level 
of many other ULXs. So, one may suspect that variability 



(which is seen in the directions of both increasing or decreas- 
ing luminosity; Fabbiano et al. 2003) occasionally transform 
ULXs in HLXs and viceversa. The dimming of source N.10 is 
such that its flux is now consistent with the star-formation- 
rate normalized XLF of the Cartwheel (assumed to be con- 
stant as observed in other well monitored sources such as 
the Antennae, Zezas et al. 2004), demonstrating continuity, 
at the brightest fluxes, with the HMXBs and ULXs. Most 
probably, then, ULXs and HLXs are low and high states 
respectively of the same class of sources. Continuity in the 
XLF with HMXB suggests also that ULXs might not be 
an entirely different phenomenon: they may represent ex- 
treme high luminous states related to the fainter accreting 
binaries. At present we can not exclude that source N.10 is 
powered by an intermediate mass black hole since we lack 
detailed information on the spectrum and on the variability 
properties due to the limited statistics available. Deep op- 
tical and radio observations might give further insight into 
the nature and the environment of such source. On the other 
hand, source N. 10 may be a HMXB accreting anisotropi- 
cally onto a stellar mass black hole in a peculiarly high lu- 
minosity state (King 2002, 2004). We can reject instead the 
possibility that a rotation powered Crab-like pulsar (Perna 
& Stella 2004) is hosted in the HLX, since a luminosity de- 
cay by magnetic braking over a time scale of six months 
would require the occurrence of a young pulsar of compara- 
ble age. This is inconsistent with the stability of the light 
curve observed over the last 4 or 10 years, from the ROSAT, 
Chandra and XMM-Newton composite data. 

King & Dehnen (2005) suggested recently that HLXs 
differ substantially from ULXs, and are naked, tidally 
stripped nuclei of dwarf galaxies hosting a massive relatively 
bright black hole. For substantial tidal stripping to occur, a 
very close distance of approach (~ 200 pc) is required for the 
impinging dwarf inside the main galaxy. This is just the dis- 
tance of X-l from the core of M82, which the authors inter- 
pret as an active relic of a naked dwarf core. For source N.10 
this scenario runs into serious difficulties. The southern ring 
of the Cartwheel is a coherent expanding wave of star forma- 
tion associated to a strong gaseous density wave excited by 
a collisional perturbation with a nearby galaxy. Given the 
gas-dynamical origin of the ring it is unrealistic to believe 
in a chance coincidence of source N.10 with a stripped core 
of a dwarf passing by. We have also investigated the possi- 
bility that we are observing a background AGN. From the 
extragalactic LF (Hasinger et al. 1993) we derive a chance 
coincidence for a background source of < 2 x 10~ 3 for the 
whole ring at the flux of the Chandra detection of source 
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N.10. This is a small but non negligible possibility, which 
however we consider unlikely. 

Individual most luminous ULXs, termed here HLXs, 
may represent occurrences of accretion episodes onto inter- 
mediate mass black holes. Source N. 10 in the Cartwheel, 
together with source X-l in M82, could be the cleanest ex- 
ample. Whereas source X-l in M82 stayed in the HLX state 
for timescales of hours, source N.10 in the Cartwheel is the 
longest lived HLX, having been observed in a time frame of 
at least 4 years to be as bright as Lo.s-iokeV ~ 1-3 X 10 41 
ergs -1 . We thus expect that other bright ULXs may share 
variability of this kind and become so bright to be ob- 
served at the HLX level. We cannot establish at the mo- 
ment whether HLXs represent an altogether physically dis- 
tinct class of objects as suggested e.g. by King & Dehnen 
(2005), or an higher luminosity transient-state of ULXs. 

Repeated observations of the Cartwheel at high resolu- 
tion are the best way to properly determine the variability 
pattern of source N.10 and the best mean to properly study 
this extreme source, which could be the first example of a 
long lived HLX in our local universe. 
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